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Feedback  systems  with  right  half-plane  poles  and  zeros  may  have  inherently 
very  poor  sensitivity  properties.  In  the  design  procedure  presented,  the 
closed-loop  poles  are  restricted  to  two  possible  regions  in  the  complex  plane. 
One  region  is  s s -a  , a>  0 . A second  is  the  interior  and  boundary  of 

a circle  in  the  left  half-plane.  The  design  is  optimum  in  the  sense  of  maxi- 
mizing the  gain  factor  uncertainty,  for  which  the  restriction  is  satisfied. 

The  design  procedure  is  very  simple  to  execute  and  results  in  loop  transmission 
poles  and  zeros  which  are  symmetrical  with  respect  to  the  boundary  of  the 
forbidden  region.  The  closed-loop  poles  lie  entirely  on  the  boundary,  over 
the  range  of  gain  uncertainty . ^ 
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1.  INTRODUCTION 


Consider  the  constrained  part  (plant)  of  a linear  time-invariant  feed- 


back system,  which  because  of  parameter  uncertainty,  has  a transfer  function 
known  only  to  be  a member  of  a set  Jp  = (P(s)}  . For  a large  class  of 
minimum-phase  (no  zeros  in  the  interior  of  the  right  half  complex  plane) 

$ sets,  any  narrow  but  nonzero  frequency-domain  performance  tolerances 
can  be  theoretically  achieved  in  a two-degree-of-f reedom  feedback  structure 
(Horowitz  1963)  even  for  very  large  but  bounded  parameter  uncertainty. 
jp  may  contain  elements  with  uncertain  right  half-plane  (denoted  by  rhp) 
poles.  A synthesis  procedure  exists,  permitting  optimum  design  to  specific- 
ations (Horowitz  and  Sidi  1972).  However  if  jP  includes  nonminimum-phase 
elements  (i.e.  with  rhp  zeros),  then  a given  set  of  performance  tolerances 
may  not  be  theoretically  achievable  V P 6 (P . It  has  been  shown  how  to 
check  if  a given  specification  set  is  achievable.  Also,  for  a given  non- 
minimum-phase but  stable  Q set,  the  problem  may  be  made  solvable  if  it  is 
permitted  to  decrease  sufficiently  the  system  bandwidth  (Horowitz  and  Sidi 
1978).  Of  course,  the  latter  solution  may  be  extremely  undesirable,  but  is 
unavoidable  in  the  linear  time  invariant  framework.  For  example,  if  3 
a P € p with  a zero  at  0.1,  the  closed-loop  bandwidth  may  have  to  be  a 
small  fraction  of  0.1  rps  i f jP  is  a 'large'  set,  in  order  to  achieve 
reasonable  tolerances  in  the  resulting  very  small  system  bandwidth.  If  ^ 
includes  both  nonminimum-phase  (denoted  by  nmp)  and  unstable  elements,  then 


even  the  latter  is  in  general  unachievable.  The  sensitivity  reduction  ■ ■ 

capabilities  of  the  feedback  loop  are  severely  restricted,  no  matter  how  White  Section 

Buff  Section 

small  the  bandwidth.  5 

1 ______ 

This  paper  considers  the  nmp  unstable  jp  problem  from  a somewhat 

different  viewpoint.  What  is  the  maximum  tolerable  plant  gain  (k)  AVAlUfiHJTT  CflOo 
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uncertainty  p = kmax/kfn . r , such  that  the  closed-loop  poles  are  restricted 
to  a certain  region  in  the  complex  plane?  One  such  region  considered  is 
^e  s s -o  , p>0  . Another  is  the  interior  and  boundary  of  a circle 
centered  at  -b  with  radius  a , b>a>0  . What  are  the  trade-offs  between 
a , b , a and  p ? A very  simple  synthesis  procedure  is  presented  for 
implementing  the  above.  All  transfer  functions  are  assumed  finite  rational 
functions.  Note  that  pure  time-delays  can  be  so  approximated  as  accurately 
as  desired,  over  any  finite  frequency  range.  Since  pole-zero  uncertainty  is 
not  considered,  it  is  assumed  that  there  are  no  hidden  cancellations  of  rhp 
poles  and  zeros  in  any  P 6 $ . 


1 . 1 One  and  two-degree-of-f reedom  structures 

A typical  tdf  (two-degree-of-f reedom)  structure  is  shown  in  Fig.  1.  The 
system  transfer  function  T(s)  = F ^ , L ■ PG  and  the  sensitivity  function 
Sp  = |yy-p-  = , are  independently  realizable  to  a large  extent  (Horowitz  and 

Sidi  1972),  because  two  independent  functions  F , G in  L = GP  are  avail- 
able to  the  designer.  In  the  odf  (one  degree  of  freedom)  system,  F=  1 , so 
T = (l-S)  . In  the  minimum-phase  ^denoted  by  mp)  system,  |s|  can  at  least 
be  designed  as  small  as  desired  over  any  finite  10  range,  making  T=1  in 
this  range,  i.e.  forcing  T to  have  a large  bandwidth.  !n  fact,  it  has  been 
shown  that  in  any  practical  design  |S|  >1  in  an  important  frequency  range 
of  T(joj)  (Horowitz  19&3)  • 

It  will  be  seen  that  in  nmp  unstable  systems,  |S|  cannot  in  general  be 
made  small  over  any  desired  a)  range.  In  fact,  |S|  tends  to  be  embarrasingly 
large  in  the  range  of  the  system  bandwidth.  In  the  tdf  structure,  one  can  at 

least  use  F to  achieve  reasonable  nominal  |T(jw)|  over  any  desired  frequency 

range,  i.e.  at  a specific  nominal  set  of  plant  parameter  values.  However,  in 
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the  odf  system  where  T = (l-$)  , not  only  does  one  have  large  |S]  , but  also 
large  | T j V P £ . The  design  challenge  is  to  do  the  best  possible  for  S . 

In  the  odf  structure  the  design  is  then  complete,  while  in  the  tdf  there  is 
available  F to  achieve  the  desired  nomimal  T(s)  . Obviously,  the  tdf 
structure  is  preferable  and  only  requires  that  any  two  independent  functions 
of  command  input  r and  output  c in  Fig.  1 be  measurable  (Horowitz  1963) • 

1.1.1  Constraints  on  S(s)  and  T(s) 

Since  rhp  poles  and  zeros  should  not  be  cancelled  for  well-known  reasons, 
the  rhp  zeros  of  P must  appear  as  zeros  of  T . Most  important,  the  rhp 
poles  and  zeros  of  P must  appear  as  such  in  L(s)  . Only  these  constraints 
need  to  be  explicitly  recognized  in  our  approach.  The  other  constraints  (rhp 
poles  of  P as  zeros  of  S , etc.)  are  automatically  included. 

1 .2  Limitations  due  to  rhp  poles  and  zeros 

Qualitatively,  the  limitations  on  the  feedback  capabilities  of  such 
systems  can  be  explained  as  follows.  If  P ( s ) has  a zero  at  b>0  , then 
the  crossover  frequency  wc  (defined  by  | L ( Jo»  ) | =1  ) in  any  practical  design 
must  be  < ab  , a<  1 (Horowitz  1963)  — assuming  that  |L(jw)|  > 1 for 
oj  < o)c  . The  latter  is,  of  course,  essential  to  achieve  the  benefits  of  feed- 
back in  a range  w < o>c  . On  the  other  hand,  if  P has  a pole  at  m>0  , it 
is  necessary  that  u>c  > Bm  , B>1  . How  can  these  two  constraints  be  satisfied 
if  both  exist  in  P with  m>b  , and  especially  if  there  are  several  rhp  poles 
and  zeros?  It  will  be  seen  that  a stable  design  with  closed  loop  poles  in  the 
restricted  regions  of  Sec.  1.1  always  exists  over  a finite  range  of  gain  un- 
certainty. The  system  overcomes  the  problem  by  having  several  crossover 
frequencies,  at  least  as  many  as  rhp  poles.  The  need  for  this  can  be  seen  also 
from  the  Nyquist  criterion,  which  requires  as  many  negative  encirclements  of  -1 
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as rhp  poles.  The  locus  of  L ( j o>)  then  tends  to  be  tightly  wound  around  -1, 
with  resulting  very  small  | 1 + L ( ju>) [ and  very  large  sensitivity  function 
S(jw)  = [1+L(jw)]  * , over  wide  w intervals. 

Consider  a typical  mp  stable  practical  L(s)  -*■  ks  x with  x *2  as 
s -*■  00  , and  | L ( j tu)  ! > 1 in  [0,uc)  . The  Nyquist  mapping  of  L(s)  (with  s 
encircling  the  rhp)  must  encircle  -m  (m  > 0)  positively,  i.e.  in  same  sense 
as  the  s encirclement,  usually  taken  as  clockwise  with  u>  : 0 -*•  <*>  . In  fact, 
it  is  impossible  for  such  L(s)  to  have  a negative  encirclement,  as  seen  from 


the  familiar  Nyquist  formula:  Number  of  encirclements  N = N -N  . Since  N 

z p p 

the  number  of  rhp  poles  of  is  0,  N = N^  the  number  of  rhp  zeros  of  L , which 


must  be  a non-negative  number.  But  when  Np > 0 , then  negative  encirclements 
are  mandatory  and  are  achievable  because  the  rhp  poles  provide  the  means,  by 


their  contribution  of  effective  phase  lead  accompanied  by  amplitude  decrease. 


But  one  must  wait,  as  uj  goes  from  0 to  «,  for  the  rhp  pole  corner  frequency 
before  the  combination  of  phase  and  amplitude  is  available  for  a negative 
encirclement.  That  is  why  in  open-loop  unstable  systems,  the  crossover 
frequency  w must  exceed  some  minimum  value  (Horowitz  1963)-  In  the  absence 
of  rhp  zeros  of  L , the  negative  encirclement  can  be  done  in  any  oi  range 


exceeding  the  minimum,  so  any  desired  benefits  of  feedback  are  achievable 
(Horowitz  and  Sidi  1972),  as  with  mp  stable  plants.  However,  a rhp  zero  corner 
frequency,  by  its  combination  of  amplitude  increase  and  effective  phase  lag, 
forces  the  completion  of  the  negative  encirclement  before  the  rhp  zero  effect 
is  too  strong  — which  is  associated  with  a maximum  value  for  u>c  . Thus,  in 
a general  sense,  lhp  poles  and  zeros  tend  to  give  positive  encirclements  while 
rhp  poles  and  zeros  give  negative  encirclements.  However,  lhp  poles  and  zeros 
can  be  effectively  cancelled  and  thus  easily  shifted,  postponing  the  encircle- 
ment (where  |l|  <1  ) to  as  high  an  u range  as  desired.  Such  cancellation 
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is  impossible  for  rhp  poles  and  zeros.  The  -1  encirclements  and  consequent 
small  |S(  = |1+L|  1 must  be  taken  in  the  range  of  occurrence  of  the  rhp  poles 
and  zeros.  In  fact,  as  will  be  seen,  it  may  even  be  necessary  to  add  more  of 
them. 

This  discussion  is,  of  course,  highly  qualitative  and  it  is  challenging  to 
make  it  quantitative.  The  technique  presented  here  does  so  to  a certain  extent. 
It  is  noted  that  contributions  to  this  problem  have  also  been  made  by  (Chang 
1961,  Brasch  and  Pearson  1970,  Shaw  1971.  Borigiorno  and  Yonla  1977).  All  of 
these  neglect  explicit  plant  uncertainty. 


2.  SOME  RESULTS  FROM  THE  BLENDING  PROBLEM 


This  paper  uses  some  results  from  the  following  "blending"  problem 
(Horowitz  and  Gera  1978):  There  are  given  two  nmp  uncertain  plants  P1  £ , 

P2  ^ $ 2 * wit^  the  same  'nPut  but  whose  outputs  can  be  separately  measured  and 
processed  by  G1  , G ^ respectively.  Can  fixed  rational  function  G1  (s)  , 

62(5)  be  found  such  that  the  sum  P ^ G ^ ^2^2  's  over  x 2 ^ 

p,g2 

Equivalently,  find  H such  that  1 + -.-r-  = 1 + PH  (with  P * P»/P,  , H = G_/G,  ) 

r j ^ L I L 1 


ha 


s no  rhp  zeros  for  any  P € $ 1 *^2  * 


If  there  is  uncertainty  only  in  the  gain  factor  k of  P (none  in  its 

Doles  and  zeros),  then  the  left  half-plane  (denoted  by  Ihp)  poles  and  zeros  of 

P can  be  cancelled  out  by  H . Those  in  the  rhp  should  not,  of  course,  be  so 

eliminated.  The  optimum  PH  (defined  as  that  which  maximizes  p“k  /k  . ) 

max  m i n 

has  the  form  PH  - kKtp(s)ip(-s)  , <p(s)  monic,  so  that  the  poles  and  zeros  of 
PH  are  symmetrical  with  respect  to  the  ju  axis.  Also,  PH  has  an  equal 
number  of  poles  and  zeros,  which  is  acceptable  because  H - G^/G^  . 


For 
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k £ [k  . ,k  ] the  root  loci  of  1 + PH  = 1 + kcp(s)ip(-s)  , all  lie  on  the  jw 
min  max 

axis  and  cover  it  completely  (see  Fig.  2).  There  is  a pair  of  zeros  of 

1 + PH  at  the  origin  either  at  k . or  at  k and  similarly  a pair  at 

min  max 

tj00  . Both  at  k . and  at  k the  other  zeros  of  (l  + PH)  are  in 
J min  max 

2 2 2 

coincident  pairs  on  the  jio  axis,  i.e.  have  factors  (s  + io.)  , as  in 

Fig.  2.  A method  was  given  for  finding  the  optimum  PH  (which  maximizes 

k „ /k  . = p ),  but  a much  better  method  is  now  presented.  The  directions 

max  min 

are  first  given  together  with  a numerical  example,  followed  by  its  derivation. 

1 

2 . 1 Directions  for  finding  optimum  PH 

Step  1 . Only  the  rhp  poles  ( 6p  in  number)  and  zeros  ( 6^  in  number) 

of  P are  explicitly  displayed  here.  If  the  system  is  to  be  type  m,  let  P 

have  m poles  at  the  origin  included  in  the  above.  The  optimum  H has 
26^  + 6z  - 2 zeros  and  6p  + 26z  - 2 poles,  (not  counting  those  which  cancel  the 
lhp  poles  and  zeros  of  P ).  Hence,  PH  has  a total  2(6p  + 6z*1)  of  poles  and 
of  zeros,  symmetrical  with  respect  to  the  ja>  axis,  viz  PH  = kK(p(s)tp(-s)  . 

Accordingly,  choose  ip(s)  with  only  lhp  poles  and  zeros  of  the  required 
generality,  i.e.  with  6p  - 1 free  (unspecified  as  yet)  zeros,  6^  - 1 free 
poles  and  the  constrained  mirror  image  of  the  6p  rhp  poles  and  the  6? 
zeros  of  P . 

Example.  P = — S-  has  6=2  . 6=4  so  let 

(s-2) (s-3)  (s  -s+4)  z p 

<p(5) » 44 , (o 

d(s)  (s+2) (s+3) (s2+s+4) (s+a) 

Step  2.  Expand  n(s)  and  d(s)  of  (1)  as  follows,  with  e. (s)  , o;  (s) 

even  and  odd  monic  polynomials,  respectively. 
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n(s)  = °n(s)  + cnen(s) 

+ (A+5) [sk  + s2 

d (s)  = od(s)  + cded(s) 


[s5  + s3(B  + 5A  + 4)  + s (5C  + 4B)  ] 

(C  + 5B  + 4A)  _4C_, 

A+5  A+5 

[s5  + s3(6a  + 15)  + s(26a  + 24)] 


(2a) 


+ (a*6)[sk  + 


.2  (15a + 26)  . 249l 
5 "Ta+T)  " + 


(2b) 


If  + 6p  ~ 1 is  even,  the  expressions  each  have  the  form  e. (s)  + c.o. (s) 
with  degree  of  e.  = (degree  of  o.)  + 1 . 


Step  3-  Make  °n(s)  = od(s) 


en(s)  = ed(s)  by  equating  coefficients, 


givin  precisely  6z  + 6p  - 2 equations  in  as  many  unknowns.  Thus,  set 


B + 5A  + 4 = 6a +15 
5a  + 

A+5  a+6 


5C  + 4B  = 26a  + 24  , 

4C  = 2 4a 
A+5  a+6 


C + 5B  + 4A  15a + 26  4C  24 a ...  , 

. — 3 — . ^ — - = — 7-  with  solutions 


a = .66006  , C = 4.14477  , A = 1.97019  , B = 5-10942  , 

s3  + As2  + Bs  + C = (s  + 1.00133)  (s2  + .96886s + 4.13928) 


Step  4.  p 


max 


m i n 


■ (^)2= 


(?) 


of  (2a, b) 


1.0953 


Step  5-  From  (2a,b)  , Step  4 gave  ^ in  PH  = kK  j s ’ Since  P 
is  known,  H can  be  found.  The  Ihp  poles  and  zeros  of  P , if  any,  must  be 
added  as  zeros  and  poles  of  H , respectively.  The  root  loci  of  1+PH  are 
shown  in  Fig.  2. 


Step  6.  The  zeros  of  1+PH  at  k . and  k are  available,  if  desired, 
c — min  max 

from  o.  , e.  of  Step  2.  Thus,  o(s)  - sU*4  + 18. 96036s2  + 41.16155) 

= s(s2  + 2.50076) (s2  + 16.45960)  , so  at  k . 1+PH  has  zeros  at 

m 1 n 

+j/2. 50076  , ±j/l6.4596  . Since  e(s)  *=  (s2  + 4.90561)  (s2  + .48487)  , at 


i 


I 


( 
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k 1+PH  has  zeros  at  ±j /A .90561  , ±j /. 43487  . 

max  J 

Derivation 

The  derivation  of  the  above  procedure  starts  with  the  results  from 

(Horowitz  and  Gera  1978)  stated  in  second  paragraph  of  Sec.  2,  giving 

PH  = kK  and  Num.  ( 1 +PH)  = d(s)d(-s)  + kKn(s)n(-s)  = 0 Kk ) -rr  ( s 2 + cj  ? ) 

2 

with  b>.  >.  0 for  k € [k  . ,k  ] . However,  at  k . , say 

j min’  max  ’ min 

d(s)s(-s)  + k^.^  Kn(s)n(-s)  = Kj  s2tt  (s2  + ok)  2 (3a) 

wh i le  at  k . say 
max  ' 

d(s)d(-s)  + kfflaxKn(s)n(-s)  = K2ir(s2+fi2)2  (3b) 

Actually  s*  may  appear  as  a factor  in  (3b)  rather  than  in  (3a)  but  it  is  not 

important  which.  K,  , K.  are  used  in  (3a, b)  instead  of  1+k  . K , 1+k  K 
12  min  max 

because  one  of  the  latter  must  be  zero,  at  which  point  a pair  of  zeros  of 

1+PH  is  at  ±j“  . Eliminating  in  turn  d(s)d(-s)  and  n(s)n(-s)  from  (3a, b) 

2 2 

give  the  pair  of  equations  (with  fit  , w.  >0  ). 

(k  " k . )Kn(s)n(-s)  = K w(s2+£3?)  - K.  s2tt(s2  + uk  ) 2 (4a) 

max  mm  / 1 l 1 

(k  . - k )Kd(s)d(-s)  = k . K,rr(s2+n?)2  - k K. s2tt (s2  + oj2)  (4b) 

' min  max'  v ' v ' min  2 1 max  1 1 


n(s)  = it(s  + z.)  , d(s)  = ir(s  + p.) 


ith  ^ e(z.)  , $e(p.)  > 0 . 


(5a, b) 


Replace  sZ  by  w so  that  (4a, b)  become 


2 2 2 2 2 2 
Y^Tr(w-Zj)  = ir(w  + fij)  “ p^W7r(w  + Wj) 

Y2it(w-o2)  = ir(w  + fl2)2  - p2wtt(w  + w2)2 


(k  - k . )K 
max  min 


(k  . - k ) K 
mm  max 

k . K_ 
mm  2 


(6a-d) 
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p 


2 

1 


2 

°2 


k K, 

max  1 

k . K0 
mm  2 


P 


O 


k 

max 

k . 
mm 


(6e-g) 


Each  of  (6a, b)  has  precisely  the  form  of  the  optimum  solution  for  active 
RC  synthesis  by  means  of  "negative  impedance  conversion  (NIC)"  (Horowitz  1959) 
There  the  problem  is  to  write  any  given  real  coefficient  polynomial  ij/(w)  , as 
the  difference  of  two  polynomials  A(w)  , B(w)  where  zeros  are  all  negative 
real,  i(i(w)  = A(w)  - aB(w)  . The  minus  a can  be  implemented  by  means  of  a 
"negative  impedance  converter".  Of  the  infinitude  of  A(w)  , B(w)  available 
for  the  task,  it  is  desired  to  choose  that  which  minimized  the  sensitivity  of 
ijj(jio)  to  variations  in  a , S^(jw)  = . It  was  shown  that  the  unique 

optimum  A(w)  , B(w)  have  precisely  the  form  on  the  right  of  (6a, b) . Thus 
n(s)n(-s)  and  d(s)d(-s)  with  w=s  have  identical  "optimum  NIC 
decomposition"  polynomials,  differing  only  in  the  gain  factors. 


(Calahan  I960)  found  a very  elegant  technique  for  deriving  the  optimum 
NIC  decomposition  polynomials  A(w)  , B(w)  from  a given  iji(w)  • Start  with 
ip (w)  = it (w  - r f ) = A(w)  -aB(w)  andwrite  N(u)*7r(u  + r.)  with  -r . the 
left  half-plane  root  of  r2  . Expand  N(u)  into  even  and  odd  polynomials. 

N(u)  * (u)  + oN(u)  = ir(u2+a.)  + puiT(u2  + b.)  , a.,b.  >0  , (7a) 

because  N(u)  is  Hurwitz  (Weinberg  1962).  Finally,  recover  iji(u)  by 
wr i t i ng 

<J>(u)  = (e2(u)  - o2(u)]  - = !T(w+aj)2  - p2w(w+b.)2  (7b) 

u^=w  1 

with  the  right  side  the  desired  optimum  NIC  decomposition. 

2 

It  follows  from  the  above  that  if  i|<1  (w)  - ir(w-zf)  of  (5a),  ^(w)  = 

2 

ir(w-p.)  of  (5b),  have  the  identical  optimum  NIC  decomposition  monic  poly- 
nomials, their  corresponding  N^  (u)  , N2(u)  in  the  Calahan  technique  also 


i 


have  identical  even  and  odd  monic  decomposition  polynomials.  The  n(s)  , d(s) 

of  (5a, b)  and  Steps  1,2  of  the  "Directions"  (Eqs.  1,2)  are  precisely  the  N. 

in  the  above  oT  n(s)n(-s)  , d(s)d(-s)  . Hence,  their  respective  monic  even 

and  odd  decompositions  must  be  identical.  In  the  Calahan  technique,  the  p 

cf  N of  (7a)  is  squared  to  recover  ip  in  (7b).  Hence  from  (6g)  , 

p = k /k  . is  obtained  by  means  of  Step  k of  the  "Directions".  Note  also 

m3  x min 

that  the  even-odd  polynomial  decomposition  of  n(s)  (or  of  d(s)  ),  gives  the 
roots  of  1+PH  at  k . , k . The  even  polynomial  gives  one  and  the  odd 

gives  the  other,  as  illustrated  in  Step  6 above. 


3.  IMPLEMENTATION  OF  OPT  I HUH  BLENDING  SOLUTION 

In  the  blending  problem,  the  boundary  of  the  undesirable  region  is  the 
jid  axis.  In  the  control  problem,  it  is  realistic  to  use  as  boundary  the 
vertical  line  s = -o  , ct> 0 , especially  as  this  permits  exploitation  of 
the  optimum  blending  results.  Clearly,  it  is  only  necessary  to  shift  the  ju> 
axis  a to  the  left  and  shift  it  back  at  the  end.  This  is  illustrated  by 
Example  2. 

3- 1 Example  2 

p _ k(s.-1 ) anj  t^e  SyStem  is  to  be  Type  1,  so  the  pole  at  the  origin  is 

s L s 

counted  as  a constrained  right  half  plane  pole.  Replace  s by  v+o  giving 

P(v)  = /— -VT7--1'+7n"l"  i and  then  follow  the  steps  of  Sec.  2.1.  Thus, 

(v-o) [ v - (2+o) ] K 

Step  1 : n(v)  = (v+1+o) (v+z)  , d(v)  = (v+o) (v+2+ct) 

Step  2:  n(v)  = [v2  + z(l+c)]  + (l+o+z)v 
d(v)  = [v2  + o(2+o)]  + (2+2o)v 
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Steps  3, **:  z(l+o)  = o(2+a) 


k 

max 

k . 
mi  n 


( 2+2a 
^ 1+a+z 


2(1+a)2 
2p2  + 4a  + 1 


and 


PH  ( v) 


[v2  - (l+o)2J(v2  - Z2) 
(v2  - O2) [v2  - (2+o)2] 


Next  replace  v by  s+a  , giving  PH(s) 


k K ( s - 1 ) (s+2a+1  ) (s+o-z)  (s+o+z) 
s(s+2a) (s-2) (s+2+2o) 


with  root  loci  of  1+PH  qualitatively  shown  in  Pig.  3 for  the  case  o=1  , 

k 


for  wh i ch  z = 1 . 5 , ■; = p = 1 . 306 1 . At  (kK)  . = - 1 there  is  a pair  of 

’ k . min 

m i n 

roots  at  -1  and  another  at  -1  ±j®  . At  (kK)  = - 1 . 306 1 , there  is  a double 

m3  x 

pair  at  -1  ± j /J  (i.e.  at  v = ±j /a (2+a)  or  s = -1  ± j /a(2+o)  ).  Note  that 
for  k € [k  . ,k  ] , the  system  poles  lie  on  the  lire  o = -1  , but  the  system 

is  stable  for  the  larger  k range  of  1.453- 


3 . 2 Far-off  poles 

In  the  optimum  blending  solution,  as  noted  PH  is  finite  as  s -*•  <*>  , 
which  is  acceptable.  But  it  is  not  acceptable  in  the  control  problem  where 
one  should  ensure  H-<-0  as  s-*-®.  In  the  above  example  it  is  therefore 
necessary  to  add  at  least  two  poles  to  H . If  these  are  far-off,  they  will 
have  little  effect  on  the  previous  "optimum"  results.  For  example,  if  they 
are  inserted  quite  close  in  at  29,  *31  the  symmetry  is  preserved  but  p de- 
creases to  1.16A  for  the  loci  to  lie  on  the  line  s = -1  . If  the  far-off 
poles  are  inserted  at  99»  *101  p = 13-06/10.35  * 1 *26  for  the  poles  on  s = -1  , 
and  I.Al  for  stability.  The  root  loci  are  sketched  in  Fig.  4a.  The  further 
off  the  far-off  poles  the  closer  p is  to  the  supremum  of  1.306.  For  some 
range  of  k the  complex  pole  damping  factors  are  quite  small.  This  can  be 
remedied  by  using  larger  a at  the  expense  of  smaller  p . Another  method 
is  presented  in  Sec.  4. 
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The  root  loci  for  the  "optimum"  design  for  the  same  problem  by  (Bongiorno 
and  You  la  1977)  is  shown  in  Fig.  Ab.  Here  p(stability)  is  only  1 . 1 8 and 
p = 1.10  for  the  poles  to  be  left  of  the  line  s = -.4  . In  addition  their 
compensation  ->-218  as  s -*■  °°  , which  is  of  course  impractical,  with  a far-off 
pole  at  -203.  The  mandatory  addition  of  at  least  one  more  pole  would 
significantly  reduce  p unless  it  is  extremely  far-off.  Note  that  this 
design  technique  has  no  provision  for  controlling  the  system  poles  over  some 
uncertainty  range,  and  the  design  procedure  is  much  more  complicated. 


3-2  The  far-off  pole  problem 


In  Example  2 of  Sec.  3-1,  the  far-off  pole  problem  was  easily  solved,  but 
the  same  approach  may  completely  spoil  the  design  in  other  cases,  as  in 

Example  3,  where  P(s)  = . Following  Sec.  3*1 » n(v)  = v+2+o  , 


d(v)  = v+1+ct  , so 


max 


(&)’• 


Note  how  p decreases  with  0 , which 


is  the  unsurprising  result  in  general, whereby  pole  damping  and  |S|  peaking 
is  traded  against  p . G(v)  = = K‘|~^2)‘  ^or  0=1  and 

G(s)  = . Suppose  far-off  poles  are  added  as  in  Sec.  3-1.  say  at 

-P1  , -p^  , then  it  is  easily  found  that  the  closed-loop  system  is  unstable 
V p|  >2  , Pj  > 0 for  K positive  or  negative  and  of  any  magnitude.  But  it 
is  easy  to  insert  lhp  far-off  poles  and  have  a stable  system  over  a significant 
k range,  e,g,  one  of  I . 85  for  a pair  of  poles  at  -30  and  larger  range  if  they 
are  further  off.  On  the  other  hand,  the  latter  approach  may  spoil  the  design 
as  i n 


Example  4,  where  P(s)  = . Following  Sec.  3*  I » G(s)  = — ^"^y  if  a = I . 

But  the  addition  of  even  one  lhp  pole,  no  matter  how  far-off,  gives  an  unstable 
system.  However  a pair  of  far-off  poles  at  p1  , -P2  may  be  inserted  as  in 
Example  2,  giving  a stable  design. 
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The  above  examples  reveal  the  importance  of  understand! ng  the  phenomenon 
and  of  finding  a simple  technique  to  determine  how  to  proceed  with  the 
additional  poles  needed  to  ensure  a practical  design.  To  understand  it,  simply 
make  an  approximate  qualitative  Nyquist  sketch  of  the  'optimum'  but  impractical 
GP(joj)  as  is  done  in  Figs.  5a~c  for  Examples  2-4.  It  is  known,  of  course, 
that  the  'optimum'  design  is  stable  so  the  number  of  counterclockwise  -1 
encirclements  must  equal  the  number  of  rhp  poles  of  GP  . This  makes  it  very 
easy  to  make  the  qualitative  Nyquist  sketches,  as  one  only  needs  to  know 
whether  to  begin  (at  GP(0)  ) to  the  left  or  the  right  of  -1.  (In  Fig.  5a 
the  pole  at  the  origin  is  counted  in  the  right  half-plane  and  the  rhp  boundary 
is  indented  to  its  left  as  shown  in  the  insert  in  Fig.  5a,  for  the  sake  of 
consistency.)  In  Figs.  5a, c the  Nyquist  locus  terminates  at  the  left  of  -1, 
so  lhp  poles,  which  cause  GP  -+  0 as  <u  -*■  « and  which  contribute  phase  lag, 
must  modify  the  locus  as  shown  by  the  dashed  1 i nes-upsett i ng  the  encirclement 
count,  no  matter  how  far-off  the  poles.  In  Fig.  5b  the  dashed  modification  is 
all  right  providing  the  lhp  poles  are  sufficiently  far-off.  However,  in 
Figs.  5a, c an  additional  rhp  pole  properly  placed  gives  one  more  negative  half 
encirclement,  as  shown  by  the  dotted  lines,  so  that  the  point  GP(®)  is  at  the 
right  of  -1,  permitting  thereafter  more  far-off  lhp  poles,  if  desired. 

Thus,  the  simple  rule  is  that  the  number  of  rhp  poles  must  be  such  that 
GP(<=°)  is  at  the  right  of  -1.  Only  then  can  far-off  lhp  poles  be  properly 

added.  They  can  then  be  added  in  a manner  very  similar  to  open-loop  stable 
designs  — so  as  to  have  little  effect  on  the  gain  and  phase  margins  in  the 
crossover  regions.  The  plural  is  necessary  here  because  there  are  as  many 
crossovers  as  rhp  poles. 

If  our  'optimum'  design  places  GP(®)  to  the  left  of  -1  there  is  no 
choice  but  to  insert  an  odd  number  of  additional  rhp  poles  (or  remove  an  odd 


number  if  possible),  before  any  lhp  far-off  poles  can  be  added.  Such  additional 
rhp  poles  can  always  be  inserted  sufficiently  far-off  to  have  as-small-as- 
desired  effect  on  the  p = k /k  . of  the  original  'impractical  optimum1 
design.  If  one  wishes,  they  can  be  put  closer  in,  considered  as  part  of  P(s) 
and  an  'impractical  optimum  design'  found  for  this  new  P(s)  , whose  PG(»)  is 
now  on  the  right  of  -1,  permitting  any  finite  number  of  additional  lhp  poles, 
if  properly  placed.  However,  such  added  rhp  poles  always  decrease  the  p of 
the  previous  optimum  design,  the  amount  of  decrease  being  smaller  the  further 
off  it  is.  Note  that  in  Example  2 of  Sec.  2.1,  a symmetrical  pole  pair  was 
satisfactory.  Additional  lhp  far-off  poles  may  be  then  added,  which  is  not 
possible  without  such  or  other  preliminary  change  in  the  rhp. 


3. 3 The  problem  of  highly  underaanped  closed-loop  pojes 

In  the  design  philosophy  of  this  section,  the  optimum  design  maximizes 


p = k /k  . for  which  a 1 1 the  closed- loop  poles  lie  in  Me  s * -o  . In  this 
optimum  design,  there  is  a permissible  k value  for  a system  pole  at  each 
point  on  the  line  s * -o  from  jO  to  j»  , so  that  for  some  k range  one  or 
more  closed  loop  complex  pole  pairs  will  be  very  highly  underdamped  —at  the 
higher  natural  frequencies,  of  course.  This  may  be  intolerable.  One  remedy 
is  to  decrease  the  permissible  k range.  One  might  consider  insertion  of  a 
complex  pole  or  zero  pair  suitably  on  the  line^S“-o  , in  order  to  focibly 
curtail  the  root  loci  of  1+GP  . Thus  in  Exaniple  3 of  Sec.  3-2  with  p = 1 , 

a kK  ( v q j 

consider  use  of  PG(v)  = — r — =“ , thereby  sacrificing  "optimality" 

(vZ  - 4)(vZ  +50) 

to  some  extent.  It  is  easily  found  that  for  kK  € [22.2,29*65]  with 

p * 1 . 336  , the  roots  are  on  s * -1  from  -1  + jO  to  - 1 ± j A . 88  wi th  a 

minimum  damping  factor  C - .2  . However,  if  the  'optimum'  design  is  used, 

r,  . » .48  for  the  same  value  of  p . Of  course  in  the  "sub-optimum"  design 
min  r 
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two  relatively  close-in  additional  poles  have  been  obtained,  relieving  consider- 
ably or  completely  the  far-off  pole  problem. 

Highly  underdamped  complex  poles  in  the  p range  can  also  be  avoided  by 
using  circles,  as  in  Sec.  4,  instead  of  vertical  lines  as  the  boundary  of  the 

forbidden  region.  However,  no  matter  how  one  squirms,  the  sensitivity  function 

Sp(jw)  = 0+L)  1 terlds  to  be  large  in  the  control  bandwidth  range,  because  of 
the  need,  noted  in  Sec.  1.2,  of  L(j<o)  to  negatively  encircle  -1  a sufficient 
number  of  times. 

3-4  Elements  in  the  Feedback  Return  Path 

In  Fig.  1,  the  return  path  transmission  from  C to  U is  -1,  implying 
the  sensor  has  infinite  bandwidth  with  value  1.  Its  actual  transfer  function 
M(s)  is  easily  accommodated  in  the  design  procedure  which  emerges  with  F(s)  , 
G(s)  in  order  to  realize  a T(s)  , S(s)  pair.  Suppose  M(s)  is  mp  and 
stable.  Then  the  same  T , S are  achieved  by  using  F*  , G*  with 
T = TT^P  " TTpfe  • S_1  " <1+GP)  * (I  +C*W0  . Hence,  set  G*  - G/M  , 

F*  * FG/G*  * FM  . If  M(s)  has  any  rhp  poles  and/or  zeros,  they  must  be 

considered  in  the  design  procedure  by  including  them  explicitly  in  the  loop 
transmission,  exactly  as  those  of  P(s)  . It  must  also  be  recognized  that  the 
rhp  poles  of  M must  appear  as  zeros  of  T(s)  . Hence,  if  T(s)  is  explicitly 

I 

» formulated,  these  rhp  poles  must  be  included,  as  well  as  the  rhp  zeros  of 

P(s)  . In  our  design  procedure  there  is  no  explicit  formulation  of  S(s)  . 

The  constraints  on  S(s)  due  to  the  rhp  poles  and  zeros  of  P , M are 
automatically  handled  by  the  design  procedure,  so  long  as  they  are  explicitly 
included  in  L(s)  . 
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h.  CIRCLES  AS  LOCI  OF  CLOSED-LOOP  POLES 

The  interior  of  a circle  in  the  lhp  instead  of  the  region  s <: -o  , 

may  be  used  as  the  acceptable  region  for  the  closed-loop  poles,  over  some 
finite  plant  gain  factor  range.  This  gives  an  optimum  design  in  which  all  the 
poles  lie  on  the  circle  boundary.  The  design  technique  is  just  as  simple  as 
in  Sec.  3-  The  transformation 


w = s + b ~ r m b - r - w(b+r) 

s + b + r * w - 1 


(8a ,b) 


maps  the  w-rhp  plane  into  the  interior  of  the  circle  of  radius  r , centered 
at  -b  in  the  s plane.  Hence  to  design,  one  first  maps  the  rhp  poles  and 
zeros  of  P(s)  into  the  w plane,  designs  in  the  w plane  in  the  manner  of 
the  blending  problem  of  Sec.  3.  and  then  maps  back  into  the  s plane.  This 
is  illustrated  by  several  examples. 


Example  3 (of  Sec.  3.2) 

P “ and  say  r = 1 , b ■ 2 so  pw  = ^ an<^  °^v ' ous ' V 

L = -k'K.  which  becomes  L(s)  = Aft.l5.  w i th  o * (-1T')^  = 

w 1 (w2-.25)  (s-1 ) (s+y)  5 

1.M  fot  the  roots  to  lie  on  the  circle.  Note  that  in  Sec.  3-2  p = (^-)  2 = 
2.25  >1.44  . However,  here  the  closed- loop  poles  are  well  damped  for  a larger 
range  — from  kK » 10/21  at  which  both  poles  are  on  the  negative  real  axis  to 
kK  *=  1 at  which  one  root  reaches  -«  . It  follows  from  the  discussion  in 
Sec.  3-2  that  lhp  far-off  poles  may  be  added,  with  no  rhp  modification  needed. 


Example  2 (Sec.  3*1)  with  b»2  , r ■ 1 


P 


k ( s- 1 ) 

rfet 


becomes 


write  n(w)  = (w+.5)(w+m) 
(w2  + .2)  + .933w  . Hence, 


P = — T — , |n  manner  of  Sec.  3.1 
W (w  - •-)  (w  - .6) 

2 * i 

* (w  + .5m)  + (.5  + m)w  , d(w)  - (w  + y)  (w  + .6)  = 
. 5m  * .2  and  p =■  t - 933/ ( - 5-Hn)  ] 2 * 1.075  for  the 
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roots  to  lie  on  the  circle.  The  compensation  has  zeros  at  ±.h  which 

map  in  the  s plane  (Eq.  8b)  as  zeros  of  G(s)  at  1/3,  -11/7-  The  far-off 

pole  situation  is  unaltered,  i.e.  it  is  necessary  to  add  at  least  one  more 

k'  (w  - 5)  , 

rhp  pole.  To  do  this,  let  P = ■. i so  in  the  procedure  of 

W (w  - y)  (w  - .6)  (w  - A) 

Sec.  3-1,  n(w)  = (w  + . 5)  (w2  + Mw  + N)  , d(w)  = (w  + y)  (w  + .6)  (w  + A)  . Expanding 
each  into  even  and  odd  polynomials  and  equating  coefficients  etc.,  gives 
(for  A =100  ) M=  lO^.Ul*  , N = 41.466  . Transforming  back  to  the  s plane, 
g i ves 

-kK(s-l)  (s  + |)  (s  + 2.9809)  (s  + 3-0195)  (s  + 1 .5711)  (s  - .33167) 

L(s)  = s(s+  1.5)  (s-2)(s  + 1.75)  (s  + 3-0202)  (s  + 2. 9602) 

with  the  root  loci  of  1+L^  shown  in  Fig.  6a.  The  new  value  of  p is  almost 

exactly  the  same  as  the  previous,  because  A was  taken  so  large.  It  is  now 

possible  to  add  lhp  far-off  poles  to  L(s)  . Note  the  reduction  in  p but  in 

return  the  system  poles  are  very  well  damped  over  this  range.  It  is  found  that 

| S ( jui)  | is  1 5 . 9db  (for  kK  = 1 ),  compared  to  a peak  value  of  2**.6db  in  the 

m3  x 

optimum  (Bongiorno  and  Youla  1977)  design.  Thus,  this  approach  gives  the 
designer  a flexible  means  of  trade-off  between  p and  • 


Example  5 wi th  b = 7 , r = 6 


P = ^7'S— vv  and  a more  conservative  circle  is  used  (Fig.  6b). 
s(s-1) 


become  w = 


s+1 

s+13 


- — ^-7^  so  the  zero  and  poles  of  P(s)  map 

W“  1 


1/13.  1/7  in  the  w plane.  A rhp  pole  at  A is  added  to  P(s)  , 
the  far-off  pole  problem.  Following  the  technique  of  Sec.  3-1* 


Eqs.  (8a, b) 
into  .2, 
to  handle 


n (s)  ■ (w  + .2)  (w2  + Mw  + B) 

■ w3  + w ( B + . 2M ) + (.2  + M)  (w2  + -2-^) 

d(s)  - (w  + —)  (w  + y)  (w+A) 

- w3  + w(  .01 1 + ,2192A)  + (.2198  + A)  (w2 +-$£_) 
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Let  B + .2M  -*  .011  + .22A  , , giving  (for  A = 100  ) B = 4.  7398  , 

M = 86.246  , p = = 1.344  for  the  loci  to  remain  on  the  circle. 

Mapping  back  into  the  s-plane  gives  a compensation  network  G(s)  with  zeros 
at  -12.8624,  -13.1409,  -1.6255,  -.3017,  ~3  and  poles  at  -13-1212,  -12.8812, 
-1.85714,  -2.5-  Left  half-plane  far-off  poles  may  now  be  added.  The  root  loci 
of  1+L  are  sketched  in  Fig.  6b. 


i 

l 


5.  INTERLACING  PROPERTY  OF  RHP  REAL  POLES  AND  ZEROS 


It  is  noted  that  in  all  the  examples,  the  loop  transmission  L(s)  which 

emerged  had  an  even  number  of  (zer°s)  between  any  two  (po'es)  on  the  positive 

po 1 63  Z6  ros 

real  axis.  This  necessary  and  sufficient  condition  for  system  design  was 
previously  noted  by  (Youl a et  al  1974)  . 1 1 i s not  expl  r c • t ly  needed  in  our  des ign 

approach  as  it  automatically  emerges.  Nevertheless,  it  is  worth  noting  and 
presenting  herewith  a very  simple  proof  of  its  necessity.  Its  sufficiency  is 
obvious  from  the  constructive  nature  of  the  design  technique. 


, . , KN(s)f(s) 

Let  L = -D7sTg(s)-  * 50  that 


Num.(l+L)  =*  Dg  + KNf  = ip(s)  , (9) 

w^ere  <p(s)  and  the  monic  polynomials  D , N , have  no  positive  real  zeros 
and  the  monic  polynomials  f , g have  only  positive  real  zeros.  Consider 
any  two  consecutive  zeros  of  g(s)  , say  at  p.  , p.  >0  . From  (9), 
ip(p.)  ■ KNf (p. ) and  <p(p.)  = KNf(pj)  . Since  N(s)  has  no  positive  real 
zeros  and  K is  a constant,  the  difference  in  signs  of  <p(p.)  and  tp(pj) 
is  precisely  equal  to  the  difference  in  signs  of  f(p.)  , f ( Pj ) . There  is 
an  alteration  in  the  latter  iff  f(s)  has  an  odd  number  of  zeros  between 
Pj  , Pj  . But  since  ip(s)  has  no  positive  real  zeros,  there  can  be  no 
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difference  in  signs  of  tp(  p . ) and  cp(  p j ) . Hence,  f(s)  must  have  an  even 
number  of  zeros  between  p.  , pj  . A similar  condition  for  the  zeros  of 
g(s)  is  proven  in  the  same  manner. 

Suppose  L(s)  has  an  excess  eL > 0 of  poles  over  zeros,  so  that  as 
s -*■  oo  f ip(s)  -*•  Dg(s)  >0  because  D , g and  in  this  case  <p  , are  monic 
polynomials.  Consider  the  largest  positive  real  zero  of  f(s)  say  at  z , 
so  D(z)g(z)  = tp(z)  and  the  sign  of  KNf  is  the  same  for  s real  > z . 

Let  g(s)  have  m positive  zeros  in  [ z ,°°)  with  the  largest  at  p . Hence, 
sgn  Dg(s)  (for  s real  >p  ) - (-l)m  Sgn  D(z)g(z)  = (-1)m  sgn  ip(z)  . But 
sgn  Dg(s)  = sign  cp( s ) > 0 for  s > p because  Og  -+■  as  s ■+■  <*>  and  ip(s) 
is  monic.  Hence  m must  be  even. 

6.  CONCLUSIONS 

A very  simple,  straightforward  design  procedure  has  been  presented  for 
feedback  systems  with  constrained  rhp  poles  and  zeros.  The  inherent 
limitations  in  the  sensitivity  reduction  properties,  which  are  difficult  to 
cope  with  in  the  frequency  domain,  are  handled  by  a design  which  restricts 
the  closed-loop  poles  to  certain  specific  regions  in  the  complex  plane.  The 
design  is  optimum  in  the  sense  of  giving  the  maximum  gain-factor  uncertainty 
for  which  such  restrictions  are  satisfied.  The  most  significant  shortcoming 
of  the  design  technique  is  that  uncertainty  in  plant  poles  and  zeros  is  not 
cons i dered . 
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Symbols 


1 h p 

left  half-plane 

mp 

minimum-phase 

nmp 

nonmi n i mum-phase 

odf 

one-deg  ree-of- f reedom 

rhp 

right  half-plane 

tdf 

two-deg  ree-of- f reedom 
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Fig.  1 . 


Figure  Ti ties 


Two-degree-of-f reedom  structure 


L.PG  , T(s)  ( F-l 


4 

in  one-degree-of-f reedom  structure) 


Fig.  2.  Example  1 from  Blending  Problem  — Root  Loci  of  1+L = 0 . 

I = -10  Vs2  - 1)  (s2  - 16)  (s2  - ( .0276) 2)  (s2  t 1. 81  Is  + 5-902) 
(s2  - <i)  (s2  - 9)  (s2  - (0.239) 2)  (s2  ± s + 4) 

Fig.  3-  Example  2 - Root  Loci  of  1 - k '-1  .Ul ~ • $). (s .+  2:})  ( s.+? ) 

s (s2  - A) (s+4) 


Fig.  A.  Comparisons  of  designs;  root  loci  of  1+L*0  . 

(a)  L = -kK(s  - .5)  (s-1)  (s  +2.5)  (s+3) 

s(s2  - k) (s+A) (s  - 99) (s  + 101) 

(b)  L = -k(s-l) (218s3  + 3938s2  - 2200s  - 1000) 

s (s-2) (s3  - 1 85s 2 - 3759s  - 50) 


Fig.  5-  Qualitative  Nyquist  sketches  (a)  Example  2 (b)  Ex.  3 (c)  Ex.  k. 

Dashed  lines  show  effect  of  Ihp  far-off  poles.  Dotted  lines  in 
(a,c)  show  effect  of  additional  rhp  pole. 

Fig.  6.  Root  loci  for  designs  with  roots  on  circle 

(a)  P = g'l'g-Y)'  i centre  at  -2,  radius  1. 

k(s-2)  . . j.  , 

s'('s~lj  ’ centre  at  -7,  radius  6. 


(b)  P 


■■  i — 

-! 


0 


Piq.  6 . Root  loci  for  designs  with  roots  on  circle 

(a)  P *»  k(s:.0.  1 centre  at  -2,  radius  1. 

(b)  P - ; centre  at  -7t  radius  6. 
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